Astronomy & Astrophysics manuscript no. 4260 


© ESO 2008 


February 5, 2008 





2 



Phase transitions in rotating neutron stars cores: back bending, 
stability, corequakes, and pulsar timing 

J. L. Zdunik^, M. Bejger^, P. Haensel^, and E. Gourgoulhon^ 

N. Copernicus Astronomical Center, Polish Academy of Sciences, Bartycka 18, PL-00-716 Warszawa, Poland 
LUTH, UMR 8102 du CNRS, Observatoire de Paris, F-92195 Meudon Cedex, France 
' jlz@camk.edu.pl, bejger@camk.edu.pl, haensel@caink.edu.pl, Eric.Gourgoulhon@obspm.fr 

o • 

O ' Received: 28.09.2005, Accepted: 06.01.2006 

<N : 

. ABSTRACT 

^ : 

, Aims. We analyze potentially observable phenomena during spin evolution of isolated pulsars, such as back bending and 
l/^ ' corequakes resulting from instabilities, which could result from phase transitions in neutron star cores. 

fSJ ' Methods. We study these aspects of spin evolution of isolated compact stars by means of analytical models of equations of state, 
, for both constant-pressure phase transitions and the transitions through the mixed-phase region. We use high-precision 2-D 
' multi-domain spectral code LORENE for the calculation of the evolutionary sequences of rotating neutron stars. This allows us 
^ \ to search the parameter space for possible instability regions, and possible changes in the stability character of rotating stars 
with phase transitions in their cores. 

Results. We determine the conditions on the density jump in constant-pressure phase transitions which leads to the back bending 
phenomena or to the existence of the unstable segments in the evolutionary sequences of spinning down isolated normal neutron 
' stars. We formulate the conjectures concerning the existence of two disjoint families of non-rotating and rotating stationary 
(y-^ ' configurations of neutron stars. To clarify the effect of rotation on the stability of neutron star we present the particular 
, case of EOSs leading to marginal instability of static and rotating configurations: marginal instability point in non-rotating 
configurations continues to exist in all evolutionary spin-down tracks. We discuss the fate of rotating stars entering the region of 
instability calculating the change in radius, energy release, and spin-up associated with the corequake in rotating neutron star, 
triggered by the instability. The energy release is found to be very weakly dependent on the angular momentum of collapsing star. 
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1. Introduction ply information on the properties of matter at density 

«; ■ .... ^ 10^^ g cm~^ expected at neutron star center. We can 

fTi ' The properties of the hieh-density matter in compact stars i i u j-- -c i ^ -n -j 

^ ^ ^ ^ *= only hope that observations ol neutron stars will provide 

' are now studied by means both of nuclear physics and as- •,, ^ • ^ u- u -n ui ^ i j- j. 

/ \ us with constraints which will enable us to select a correct 

trophysics. Many different equations of state (EOS) were ^^^^^ ^^^^^^ ^^^^^j^ j^^^^ ^.^^ 

proposed in the literature, based on different theoretical ^^.gptable dense matter theories, 
models of dense matter. In particular, it has been postu- 
lated that at high enough density, matter contains hy- Particularly interesting method of searching for the 
perons, and that phase transitions to exotic phases of P^asc transition in neutron star cores vianulsar timing 
dense matter, such as meson (pion, kaon) condensates, or ^^s proposed by I Cdendenning et al. | \]MW- As a pul- 
quark matter can occur in the dense compact star cores ^P"^^ down, its central density increases, and for a 
(for a recent review o f possible stat es of matter in com- certain density a new phase of matter can appear. In 
pact star cores see e.g. lGlendenni^EoOfil) . Unfortunately, ^^e case considered by | Glendenning et al. |, the new phase 
the lack of precise knowledge of strong interactions be- consisted of quark matter. The authors suggested, that 
tween hadrons in dense matter, as well as deficiencies and ^^e softening of the EOS, induced by the formation of 
approximations plaguing the many-body calculations of ^^e new dense phase, leads to a temporary spin-up era, 
the EOS prevent us from knowing the actual structure of ^^e phenomenon called back-bendmg. Originally, the name 
neutron star cores. Terrestrial experiments cannot sup- comes from nuclear physics, where the phenomenon of 

"back-bending" was observed in the systcmatics of the 

Send offprint requests to: J.L. Zdunik moment of inertia of excited states of rapidly rotating 
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nuclei, see e.g., iRin g fc Schiici] l|l98(]() . The calculations 
of iGlenden'ning et all we re performed within the slow- 
rotat ion approximation (jHartld Il967t iHartle fc Thornj 
Il968|) . supplemented with additional relations resulting 
from accounting for the rotational stretching and fr ame- 
dragging effects (IWeber fc Glendennind.ll99lLll992l). 

Se veral oth er authors |Heiselbergfc^iorth-Hens'enl 
Il998l iGhubarian et all. l200nir carried out their calcula- 
tions of the back-bending phenomenon using the slow ro- 
tation approx i matio n of Hartle. However, as shown by 
ISalgado et al.l ( 19941) . the Hartle method, when compared 
with results of exact 2-D numerical codes, breaks down for 
angular velocity close to the Kcplcrian one. 

First calculations concerning back-be nding based 
on ex act 2-D code were performed by ICheng et al.l 
1I2OO2I) . These authors used the version of KEH 
cod e (iKomatsu. E ri guchi fc HachisijL Il989aib|) . improved 
by IStergioulas fc FriedmanI lll995l) (see also references 
therein). In their work. ICheng et all focused on the role 
of the crust for the very existence of the back-bending. 
Indeed, as they show, even a slight change in the physical 
state of the crust (for example, a change in the crust-core 
transition pressure) may significantly affect the results. 
This shows that high precision is mandatory for reliable 
calculation of the back-bending phenomenon. 

Another important article c ontaining results based 
on 2- D computations with the IStergioiilas fc FriedmaE 
code, was published by ISovrou fc Stergioulas 
They showed t hat the results obtained by 
IGlendenning et al ] ^1997^ are plagued by large numeri- 
cal uncertain ties. For example, the very same EOS as 
that used by IGlendenning et all (EOS from Table 9.2 of 
iGlendennind l200(jl did not yield the back-be nding phe- 
nomenon at all in lSpvrou fc Stergioulaj l)2002(l ! It became 
evident that the back-bending problem is much subtler 
than previously considered, and that it requires careful 
handling as well as hig h-precision 2-D computations. Also, 
ISpvrou fc StergioulasI pointed out some errors in previ- 
ous papers on back-bending. For example, the formula 
for braking index must be corrected by taking into ac- 
count the rotational flatte ning of the star (Sect. 6 in 

ISwrou fc Stergioula jEmi. 

Most recently, Izdunik et aP l|2004l) showed that the 
back-bending phenomenon can also occur for the EOSs 
different from the mixed-phase one. They also pointed 
out importance of the stability with respect to the axi- 
symmetric perturbat ions. The appearance of hyperons 
in the dense matter ( Balberg fc Gal Il997l) softens some 
EOSs so much, that pulsars losing angular momentum 
actually spin up during a period of time. Paradoxically, 
during this spin-up phase pulsars could lose a significant 
amount of their angular momentum. 

It should be str essed that all previous works (except of 
IZdunik et al1l2004l) considered back-bending as a feature 
of the the dependence of the moment of inertia, /, on the 
rotation frequenc y, f = 1 /rotation p eriod: / = /(/). As 
shown recently bv lZdunik et al this can easily lead 

to incorrect determination of the stability of rotating stars. 



IZdunik et al.l l)2004|) pointed out that many cases claimed 
before to correspond to the back-bending, actually cannot 
be realized in nature because of the instability with respect 
to the axi-symmetric perturbations. 

One of the aims of the present work is to determine 
reliably and precisely the stability regions on the back 
bending segments of the spin evolution tracks. In order 
to avoid any precision problem, and to investigate large 
and possibly complete parameter space, we will work with 
analytical EOSs of dense matter exhibiting a softening 
at supra-nuclear density. Two examples of softening by a 
phase transition will be considered. We will study EOSs 
with constant-pressure phase transitions, characterized by 
a density jump obtained using the Maxwell construction. 
Second, we will use EOSs with phase transition extending 
over a finite pressure range in which two pure phases coex- 
ist forming a mixed-phase state. Such EOSs are obtained 
for the first-order phase transition between two pure phase 
by relaxing the condition of local electric charge neutrality 
and replac ing it by less string ent condition of the global 
neutrality l(Glendenninglll992i . Mixed-phase state can be 
realized provided the surface tension at the interface be- 
tween the two pure phase is not too large. 

The plan of the article is as follows. In Sect. |5] we 
present various types of analytical EOSs used in the cal- 
culation of the spin evolution tracks. Numerical methods 
used in exact 2-D calculations are briefiy presented in 
Sect. O Our numerical results are described in Sect. ^ 
We first describe the general criteria for the back-bending 
and the stability for spinning-down stars. Then, the results 
for the EOSs with a mixed-phase segment are studied in 
Sect. 14.11 and those with constant pressure phase transi- 
tion with density jump are reviewed in Sect. 14.21 In Sect. 
|S1 we describe a link between the existence of unstable 
segments in the families of static and rotating configura- 
tions of neutron stars. Change in neutron star parameters, 
accompanying transitions between two rotating configura- 
tions, triggered by instabilities of isolated rotating neutron 
stars, are studied in Sect. Modifications in the pulsar 
timing and pulsar age evaluations, due to phase transi- 
tions in spinning-down isolated neutron stars are studied 
in Sect. 13 A summary of our results and their discussion 
is presented in final Sect.jSI Formulae referring to analyt- 
ical models of EOSs with phase transitions are collected 
in the Appendix. 

2. Analytical EOS 

Studies of stability of stationary rotating configurations 
require high accuracy of the 2-D calculations. This is par- 
ticularly true for the determination of a critical configu- 
ration which separates the stable configurations from the 
unstable ones. On the other hand, in our search for criteria 
which link the properties of the phase transition in dense 
matter with its effect on the stability of stellar configura- 
tions with new-phase cores we wanted to explore possibly 
large space of the phase transition parameters. In view 
of this, we decided to use analytic forms of the EOSs. In 
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Fig. 1. Two examples of EOSs employed in our calcu- 
lations: constant-pressure phase transition (dashed line, 
EOS with Ta = Tb = 2.25, the density jump A = 1.3, 
can be found in Table and transition through a mixed- 
phase state (solid line, thin dotted lines mark m and ni\ 
EOS MM of Table [TJ. Top panel: v{p) diagrams. Bottom 
panel: density profiles in stellar core, for a M = 1.33 Mq 
neutron star models. 

this way. the EOSs were easy to handle, and we could 
easily change the phase-transition parameters in a contin- 
uous way, detecting with high precision the appearance of 
instability. Moreover, with our numerical methods of solv- 
ing the hydrostationary equilibrium problem, we had no 
unnecessary preoccupation with control of precision (see 
Sect.EJ. 

2.1. Constant pressure phase transition 

Normal phase is described by a polytrope of index Fa. 
The phase transition to a dense phase takes place at a 
specific pressure, Pab, at which normal phase of baryon 
density ua coexists with dense phase of baryon density 
riB, with baryon density jump rie — tt-a > 0. Coexistence, 
equivalent to phase equilibrium, corresponds to equality 
of baryon chemical potential of both phases. The high- 



density phase B is a polytrope of index Fb. Details of 
the construction of the analytical EOSs are given in the 
Appendix lA.ll 

2.2. Phase transition via the mixed-phase state 

Normal phase is described by a polytrope of index Fa- 
Mixed phase starts at baryon density ni, and ends at 712; 
it has a polytropic EOS with index F < Fa. This approx- 
imation is rather good for a sizable lower-density part of 
the mixed-phase segment (seelBci ger et al. 2QQ5,) . To be 
specific, pure high-density phase, which exists at n > ^2, 
was assumed to be quark matter, with a bag-model EOS. 
Details of construction of the analytical EOSs are given 
in the Appendix lA. 21 Depending on the values of param- 
eters, we obtained three possible types of the EOS with 
mixed phase region. First, we obtained EOSs which give 
exclusively stable models between the minimum and maxi- 
mum allowable mass; they are called MSt (Mixed Stable). 
Second, we obtained the EOSs which produced two dis- 
tinct stable families of neutron stars, i.e., low density fam- 
ily and high density family. They are represented by two 
segments in the mass-central density plane and are sep- 
arated by a segment of unstable configurations. Such an 
EOS is denoted by MUn (Mixed Unstable). There are 
also EOSs which yield a vanishingly short unstable seg- 
ment of static configurations, so that the configurations 
are "marginally stable there" (an inflection point). This 
class of EOSs will be denoted by MM (Mixed Marginally 
Stable). Parameters of examples of analytical EOS of each 
of these three types are given in Table and their P{p) 
plots are shown in Fig. [3 These examples of the EOSs 
ha ye rti similar to the models tabulated in Tables 9.2-9.3 
bv IClendennind (|200fll) . Our v alues of n?. ~ 3ni are sys- 
tematically lower than those of lGlendennin g (20o3) which 
are characterized by 712 ^ 4ni. However, this difference 
is not relevant for the general results of our analysis of 
stability of rotating stars. 

3. Numerical methods 

Hydrostationary stellar models have been computed 
by solving the Einstein equations for stationary axi- 
symmetric space-time. The numerical calculations were 
performed by means of the rot star code, a part of the 
LORENE public-domain object-oriented C-I--I- scientific li- 
brary based on spectral methods. For the complete set 
of equations as well as the tests and the description of 
th e numerical code we re f er the reader to t he ar ticles 
bvlBonazzola et al. I lll993j) . iBonazzola et al. I l)l998|) and 
ICourgoulhon et al. Il)l999|) For the purpose of the present 
work, we have employed two (in the case of a constant- 
pressure phase transition) and three (phase transition via 
a mixed-phase state) spectral domains to describe the stel- 
lar interior and surface. The division into domains makes 
use of the adaptive coordinates while setting the bound- 
aries between different phases. Therefore, the density field 
is smooth in each domain and the resulting accuracy is 
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Table 1. Main parameters of the EOSs with mixed-phase 
segment. Below the mixed-phase transition point ni a 
polytropic EOS with Fa is used. Mixed phase extends 
within ni < 71b < '^2j a-nd is described by a polytrope with 
adiabatic index T^^. Above the density 712 we assume pure 
quark matter with MIT bag model EOS p = ^{p—po)c^. In 
all cases the dimensionless polytropic pressure coefBcient 
K was equal 0.025 (see Appendix A for details). M^*,^*^x 
and Mjjj^* denote the maximum allowable baryon and 
gravitational mass of the non-rotating star. The EOSs are 
labeled as follows: MSt produces a stable back bending, 
MUn - an unstable one, and MM produces a marginally 
stable case (for more details see the text). 
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Fig. 2. Examples of the three EOSs with phase transition 
through the mixed-phase state, considered in the present 
paper. The parameters of the EOSs are given in Tabled 

very high - in terms of GRV2 a nd GRV3 virial error indi- 
cators (see lNozawa et ai~1ll998|) the relative error was on 
average ~ 10"''. 

4. Back-bending and stability of rotating 
configurations 

As it has has been already mentioned by us in the previous 
paper treating spec ifically the back bending in hyperon 
stars ijZdunik et al.l . |^04) , the back bending phenomenon 
itself is strictly connected with the existence of a minimum 
of Mb along the sequence of rotating configurations with 
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Fig. 3. Total angular momentum versus rotation fre- 
quency / (left panels), and moment of inertia / = J/fl ver- 
sus / (right panels), for EOSs MSt and MUn. The stability 
criterion is easily applied to left panels. It is clear that for 
the MSt EOS back bending feature is not associated with 
an instability, with all configurations being stable. On the 
contrary, the MUn EOS produces back bending with a 
large segment of unstable configurations. Simultaneously, 
the /(/) curves for both EOSs are very similar and appar- 
ently show very similar back-bending shapes. 



fixed /. More generally, the softening of the EOS leads 
to the flattening of the Af(i?oq) and M{pc) curves repre- 
senting these sequences. Here, i?oq is the circumferential 
equatorial radius. The effect is larger for larger values of 
/, and it appears above some critical value of /. The soft- 
ening of the EOS results then in a local maximum of Mb. 
In such a case, there exists a region, close to the local 
maximum mass, in which the decrease of J leads to the 
increase of /, and this is equivalent to back bending. The 
fragment of the curve, for which Af decreases as a func- 
tion of pc, does not necessarily correspond to the insta- 
bility region - the decrease at fixed / does not imply the 
decrease at fixed J. Only the latter condition is equivalent 
to the instability of rotating confi guration with respect t o 
the axi-symmetric perturbations ijFriedman et al.l Il988|) . 
Summarizing, the particular shape of the M(i?eq) curve 
for a fixed / indicates the presence of the back bending 
phenomenon, whereas the shape of M{Rcq) at a fixed to- 
tal angular momentum J tells us if the configurations are 
stable. 

In the present paper we will show examples of the EOS 
with constant-pressure (density-jump) phase transitions 
and mixed-phase transitions for which one of the two sit- 
uations occurs: 
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— 1. all configurations are stable but the back bending 
exists. 

— 2. the phase transition results in the instability region 
for rotating stars. 

To distinguish between these two cases we can look at 
the behavior either of the curves Mb(pc)j or J{^l)Mi,', the 
quantity fixed along a sequence is indicated by the lower 
index. For such curves the instability criterion directly ap- 
plies. However, is not so easy to detect the instability using 
the /(/)j\fjj plot, which may be similar to that correspond- 
ing to a fully stable sequence. An example is presented in 
Fig. 121 where we plot two functions: J{f)M^ and I{f)Mt, 
for the two models of EOS (MSt and MUn) which are 
described in detail in Sect. 14.11 Upper panels correspond 
to the stable model MSt and lower to the MUn model 
for which a region of unstable configurations exists. The 
difference between the MSt and MUn cases is clearly vis- 
ible in left panels (J(/)Mb curves), where we can easily 
recognize the instability region for MUn model, by apply- 
ing the condition (dJ/dpc)j^j^ > 0, and keeping in mind 
that pc is monotonic along this curve. However, the upper 
and lower right panels {Iif)Mi,) are quite similar, without 
any qualitative difference. Simultaneously, the back bend- 
ing phenomenon in the MUn case is very large (almost by 
200 Hz), and obviously the phase transition results in an 
instability (increase of J for increasing pc). And yet, /(pc) 
is a monotonic decreasing function all the time. 

4.1. Transition to a mixed phase 

In Figs. I4I6I we present the results obtained for three dif- 
ferent choices of the EOS's parameters describing mixed 
phase phase transitions (Table ^1. 

Model MM has been chosen to be a marginal one - 
there exists a region where the curve Afb(-Rcq) (or Aib(/5c)) 
is locally horizontal, which means the marginal stability. 
The fimction 7\/b(i?cq) for non-rotating star (/ = 0) has 
an inflection point slightly above the configuration where 
phase transition to the mixed phase occurs. At a fixed /, 
the condition for an inflection point reads, 



dpc 







dpi 



= 0. (1) 

In a special static case (/ = 0) this condition is fulfilled 
for Mb ~ 1.42 Mq and i? ~ 14 km. For this marginal EOS 
model, the back bending phenomenon exists at all rota- 
tion rates, and for baryon masses larger than 1.42 M©. 
The region for which back-bending exists is bounded by 
the points corresponding to the local maxima and minima 
of the function Mb(i?oq) at fixed rotational frequency /. In 
the marginal model for which the inflection-point condi- 
tion, Eq. J^l, is fulflUed for non-rotating stars, even an in- 
finitesimally slow rotation leads to the minimum and max- 
imum of the function Mb(i?)/ (albeit infinitesimally close 
to each other), which is equivalent to the back-bending 
(see dashed curves in Fig. ^J. However, imposing rotation 
does not result in the instability - the functions Mh{pc)j 




400 600 
f [Hz] 



1000 



Fig. 4. Upper panel: stellar baryon mass vs radius for fixed 
frequency (dashed lines) and fixed total angular momen- 
tum (solid lines), for the model MM of mixed-phase EOS. 
Bottom panel: stellar angular momentum as a function of 
the rotational frequency for fixed baryon mass (indicated 
as a label, in solar masses) for the same MM EOS. This 
EOS corresponds to the marginal case from the point of 
view of stability - the curves M(i?eq)j and J(/)Mb have 
fiat horizontal regions. The regions of back-bending are 
drawn by thick lines. 



are monotonously increasing, independently of the rota- 
tion rate (solid curves in Fig.^. Strictly speaking, for each 
J there exists an infiection point of the function M\^{pc),j 
(within the accuracy of our calculation). The size of back 
bending, defined as a frequency range for which / increases 
for decreasing J (and increasing pc), depends on the mass 
of the star (see bottom panel of Fig.^. In the model MM, 
the back bending range varies from (for A/b < 1.A2Mq) 
to ~ 180 Hz for the star with mass A/b = 1.685Afo, equal 
to the maximum mass of non-rotating configurations (i.e. 
not supra-massive ones) 

Model MSt produces a set of stable configurations in 
the region of phase transition (Fig. (S)). Here, A^b is an 
increasing function of pc for non-rotating configurations, 



6 



J. L. Zdunik et al.: Rotating NSs and phase transitions 




400 600 
f [Hz] 



800 



1000 



Fig. 5. The same as in Fig.Hbut for the MSt EOS. Phase 
transition does not result in the stabihty foss - all config- 
urations are stable. 



as well as along all rotating sequences of configurations 
at fixed J, terminating at the global maximum mass con- 
figurations. The back-bending phenomenon is in this case 
limited to the baryon masses larger than that for which 
the curve A/b(-Roq) (or equivalently Mb(/9c)) at fixed / has 
a flat, horizontal region (strictly speaking, where the Eq. 

is fulfilled). The numerical values for the MSt model 
are Afb > 1.45M0 and / > 400 Hz (see for example the 
bottom panel of Fig.|3 where the curve for A/b = 1.45Afo 
has almost vertical part). 

For the MUn model there exists a region for which 
configurations are unstable, i.e. the baryon mass is de- 
creasing function of at fixed J (Fig. ^ . In some sense 
this instability is not very strong - the difference between 
the maximum (local) and minimum mass is of the order of 
0.3%. However this feature (existence of instability region) 
is characteristic to all rotational frequencies - for all val- 
ues of angular momentum (J) fixed along the curve, the 
baryon mass has local maximum and local minimum con- 
nected by an unstable sequence of stellar configurations. 




Re, [km] 




400 600 
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1000 



Fig. 6. The same as in Fig. 01 but for MUn EOS. Upper 
panel: Phase transition results in the stability loss, i.e., 
there exists a region where Afb decreases with increas- 
ing pc at fixed J (marked by dotted lines). Lower panel: 
plots in the J — f plane, where on the dotted seg- 
ments J increases as pc increases (unstable configura- 
tions). Both features indicate instability with respect to 
the axi-symmetric perturbations. 



4.2. Constant pressure phase transition 

For non-rotating configurations, the reaction of the star 
to a constant pressure (first order) phase transition has 
been studied in de tail in the second half of 1980s (see 
IZdunik et al.lll98^ and references therein). The appear- 
ance of a new, dense phase in the center of the star re- 
sults in the change of the derivatives of the global stel- 
lar parameters with respect to Pr (see the formula B6 in 
the appendix of an article bv IZdunik et alJll987^ . Two 
important dimensionless parameters are: fractional den- 
sity jump A = pb/pa and the relativistic parameter 
XA = Pab/{pac^)- There exists a critical value of A, 
Acrit — |(1 + 2;a), such that for A > Acrit configurations 
with an infinitesimally small B-phase core are unstable 
with respect to collapse into a new configuration with 
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a large core of the dense phase. Putting it differently, a 
phase transition with A > Acrit destabilizes the star at 
central pressure Pc = Pab at which the phase transition 
occurs. It should be stressed that while A < Acrit guar- 
antees stability of small-core configurations, it does not 
assure the stability of configurations with a finite, or ~ 
in an extreme case - a large core. In such a case the in- 
stability would result from the softness of the B phase 
somewhat above pe and not directly from an over-critical 
A. In other words the compressibility of a matter leads to 
the larger mean density in the core than the value pB at 
the phase boundary. The response of the whole star to the 
appearance of the dense core built of the S-phase of the 
matter is determined by the mass and radius of this core 
(strictly speaking, this statemen t is tru e for non-rotating 
configurations, see lZdunik et all l)l987j) : for rotating ones 
also rotation rate and resulting oblateness play role). As 
a result even if A < Acrit the first order phase transition 
can lead to the unstable configurations for finite size of 
the core. As a result for the given model of the matter 
in the phases A and B there exist the maximum value 
of density jump Amax for which all configurations below 
maximum mass are stable. Of course Amax ^ Acrit and the 
difference between Amax and Acrit is larger for softer EOS 
in the phase B. 

Numerical results for a collection of sets of EOSs with 
constant pressure phase transition are collected in Table[21 
The parameters presented in this table correspond to the 
onset of back bending, i.e., the rotational frequency and 
baryon mass for which the curve Mb(i?oq)/ or /(J)A/b 
starts to have a flat region. More precisely, at these val- 
ues of frequency and mass, an inflection point appears in 
the curves under consideration. We also included param- 
eters of those EOSs for which all non-rotating stars with 
Mb < are stable. This means that for such EOSs 

the Mb(pc) curve for static configurations increases mono- 
tonically up to A^b'max- "^^^ parameter Amax gives then 
the maximum value of the density jump for a fixed set of 
other EOS parameters (adiabatic indices Fa and Fb, num- 
ber density threshold tia) for which this property of neu- 
tron stars is valid; in other words Amax corresponds to the 
"marginally stable" case. Increasing A implies increasing 
softening of the EOS by the phase transition. If A > Amax, 
the phase transition leads to the existence of an unsta- 
ble branch of the non-rotating stellar configurations. This 
unstable branch separates stable family of neutron stars 
with A-phase cores from a second family of superdense 
neutron stars with B-phase cores: these are two distinct 
neutron-star families. It should be mentioned that this 
feature (existence of the unstable region) does not depend 
on rotation - the unstable branches exist also for rotating 
configurations (strictly speaking for any value of a total 
angular momentum of the star J there exist a region with 
{dMh/dpc)j < 0). We have tested this feature (existence 
of or the lack of unstable regions) for very small depar- 
tures from marginally stable case (|Amax — A| < 0.005). 
From numerical results it follows that if A < Amax all ro- 
tating configurations are stable (before loosing stability 



at maximum mass point) and if A > Amax we have two 
branches of stable configuration for rotating stars (for any 

J). 

Picking up the onset parameters is visualized in Fig. 
[7| where we display Mb(i?oq) for one of the EOSs from 
Table El The curves are plotted for three frequencies, 
with middle one corresponding to the back bending on- 
set, / = /on. Last column of Table [5] gives the maximum 
allowable baryon mass for static configurations, M^^'^^^. 
We restrict ourselves to back bending for the normal (non 
supra-massive) stars, which appears during the spin-down 
evolution which terminates eventually by a non-rotating 
stable configuration. The dependence between the back- 
bending onset parameters - / = /on, corresponding baryon 
mass, and the intrinsic parameters of the EOS - the den- 
sity jump A, as well as the "departure" from the critical 
configuration (AA = Amax — A) is presented in Fig.|Hl The 
three families of curves visualize the data from Table [21 
(sohd lines for Fa = 2, dotted for Fa=2.25 and dashed 
for Fa=2.5). The value Amax defines the onset of back- 
bending at the limit / = 0; in this case the back-bending 
phenomenon is present for any rotational frequency. As 
it can be seen on the left panel, the onset frequency /on 
depends very weakly on the EOS in the dense core - the 
main parameter describing the reaction of the star to the 
appearance of this phase transition is the density jump. 
The right panel presents the same data not normalized 
with respect to the maximum density jump Amax - the re- 
sults can be very well approximated by the dependence 
/2„ = a AA -I- 6 (AA)3/2. These two plots can be treated as 
a slice through the parameter space to search for regions 
of the back-bending appearance - in the right panel, the 
back-bending is present, for a particular model, above its 
curve. 

5. Rotation and stability/instability of normal 
configurations 

Rotation can influence stability of a star of a given baryon 
mass. In particular, stars with Mb > M^*^*^^ cannot exist 
without rotation and collapse into black hole as soon as 
their rotation frequency falls below a certain minimum 
value. Here, however, we will restrict to a different problem 
of stability, which will concern the normal configurations 
only. 

In what follows we will use the term "stability" in a re- 
stricted sense. Namely, by stability (instability) of an equi- 
librium configuration we will mean stability (instability) 
with respect to radial perturbations in the non-rotating 
case, and with respect to axi-symmetric perturbations for 
rotating configurations. 

We studied a very large set of EOSs with phase tran- 
sitions at constant pressure, as well as those with tran- 
sition through a mixed phase state. We then produced 
static sequences and normal rotating sequences for these 
EOSs. Our calculations were very precise, because we 
used analytic forms of the EOSs. The results for both 
constant pressure phase transitions, and those proceed- 
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Fig. 7. The definition of the frequency of the onset of 
back-bending phenomenon, /on, and corresponding mass, 
M^". In the presented example /on = 449 Hz and M^"" = 
1.325Mq 
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Fig. 8. The onset frequency of the back-bending phe- 
nomenon as a function of the departure of A from the 
maximum density jump Amax (Amax — A, left panel), and 
density jump A as a function of /on (right panel) for mod- 
els presented in Tabled (solid lines for Fa = 2, dotted for 
rA=2.25 and dashed for rA=2.5). 



Table 2. Selected sets {Xa, Ta, Tb, ^a, A} of EOSs with 
a constant pressure phase transition. For all cases, A'a = 
0.025 and tt-a ~ 0.25 (see the text and Appendix IA.1|) . 



EOS 


A 


/on [Hz] 


Mr [Mq\ 






1.344 


247 


1.033 


1.256 


rA=2, 
rB=2, 

Amax= 1.364 


1.318 


368 


1.071 


1.283 


1.265 
1.212 
1.159 


482 
555 
601 


1.131 
1.193 
1.248 


1.342 
1.407 
1.477 




1.106 


624 


1.286 


1.554 




1.371 


189 


1.021 


1.475 


rA=2, 


1.318 


407 


1.087 


1.534 


rB=2.25, 


1.265 


507 


1.149 


1.597 


Amax = 1.382 


1.212 


570 


1.208 


1.667 




1.159 


620 


1.272 


1.743 




1.371 


257 


1.034 


1.712 


rA=2, 


1.318 


453 


1.110 


1.776 


rB=2.5, 


1.265 


522 


1.160 


1.844 


-^max — 1 .393 


1.218 


576 


1.215 


1.918 




1.159 


608 


1.258 


1.999 




1.376 


412 


1.307 


1.642 


rA=2.25, 
rB=2.25, 

Amax = l-432 


1.322 


537 


1.387 


1.709 


1.268 
1.215 


611 

662 


1.463 
1.536 


1.782 
1.861 


1.161 


694 


1.597 


1.947 




1.107 


709 


1.636 


2.041 




1.430 


222 


1.242 


1.803 


rA=2.25, 

rB=2.5, 

Amax = 1.444 


1.376 


449 


1.325 


1.859 


1.322 
1.268 


556 
623 


1.403 
1.477 


1.940 
2.016 


1.215 


668 


1.545 


2.099 




1.161 


696 


1.602 


2.190 




1.437 


361 


1.559 


1.967 




1.382 


532 


1.661 


2.040 


rA=2.5, 


1.329 


624 


1.754 


2.119 


rB=2.5, 


1.273 


686 


1.844 


2.205 


Amax = 1.472 


1.219 


731 


1.933 


2.297 






I ( 


9 nm 

Z.UUi 






1.109 


780 


2.076 


2.509 




1.437 


428 


1.590 


2.348 


rA=2.5, 

Tb^S, 
Amax=1.487 


1.382 


566 


1.689 


2.427 


1.329 
1.273 
1.219 


646 
701 
738 


1.871 
1.870 
1.945 


2.512 
2.603 
2.703 




1.164 


761 


2.011 


2.811 



ing through mixed phase, turned out to be qualitatively 
the same. In all cases, if non-rotating configuration were 
stable (monotonically increasing Mb(pc) and M{pc)), then 
for any value of the total angular momentum J the func- 
tions Mb(pc)j and M{pc)j were monotonically increas- 
ing, too. Thus, when all non-rotating configuration with 
Mb < M^^^^^ were stable (with respect to radial pertur- 
bations), all normal rotating configurations were stable 
too (with respect to axi-symmetric perturbations) . On the 
other hand, if for non-rotating stars there existed a region 
with decreasing Mb(pc) and M{pc), even extremely small 
one with a very shallow minimum, then an unstable region 



persisted within the rotating configurations, at each value 
of J. These two cases are illustrated in the Figs. |5l and El 

We studied also the case of marginally stable EOS. An 
inflection point, witnessing marginal stability, present in 
the Mb(i?) or M{pc) curves for non-rotating stars contin- 
ued to exist in the Mb(i?oq) or M{pc)j curves for normal 
rotating stars (Figs^EJ). 

The analysis of numerical results leads us to an in- 
teresting conclusion. Namely, for an EOS with a phase 
transition (constant pressure one or through mixed-phase 
state), rotation neither stabilizes nor destabilizes normal 
sequences of stationary configurations based on this EOS. 
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Fig. 9. The mass M of the star as a function of central 
baryon number density ric for the MM model. Dashed hues 
- fixed rotational frequency. Solid lines - fixed total angular 
momentum J . The vertical lines correspond to the densi- 
ties of transitions from the nuclear to mixed phase and 
from the mixed phase to the pure denser phase (here - 
quark matter) - see Table 

We define a family of configurations as a compact set 
(in mathematical sense) of configurations. Similarly, an 
EOS leading to a marginally stable point for non-rotating 
stars, produces also spin-down evolution tracks with a 
marginally stable point. Our result can be formulated as 
three conjectures: 

I J/ an EOS with a phase transition gives a single family 
of static stable neutron stars then it produces also a single 
family of stable rigidly rotating normal stars. 

II // an EOS with a phase transition gives two disjoint 
families of stable static stars then it gives also two disjoint 
families of stable rigidly rotating normal stars. 

III // an EOS with a phase transition gives two fami- 
lies of stable static stars separated by a marginally sta- 
ble configuration, then it gives also two families of stable 
rigidly rotating normal stars separated by a line consisting 
of marginally stable configurations. 

6. Corequake resulting from instability 

As was shown on many occasions in the preceding sections, 
an isolated neutron star, loosing its angular momentum, 
moves down along the line of fixed Mb in the J(/) plane, 
and can at some moment reach the instability point (i) ( 
the minimum of J(/) at fixed Mb in Fig. If The sub- 
sequent behavior of the star cannot be described by our 
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Fig. 10. Evolution of an isolated pulsar loosing angular 
momentum, after it reaches the instability region in J — / 
plane and then collapses. Arrows lead from unstable con- 
figuration to a collapsed stable one, with the same baryon 
mass and angular momentum. Dotted lines - unstable con- 
figurations. 
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Fig. 11. Changes of stellar parameters of a rotating soli- 
tary neutron star, due to a collapse which occurs after 
a pulsar loosing angular momentum reaches an unstable 
configuration. 
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Fig. 12. Total gravitational mass of the star as a func- 
tion of its angular momentum, for fixed baryon number 
of the star for our EOS model MUn. The central density 
is increasing along this curve as marked by the arrows. 
The upper segment (dotted) corresponds to the unstable 
configurations. Two cusps reflect strict property that the 
mass and angular momentum have simultaneous extrema 
along the path with fixed baryon number. 



stationary rigidly rotating model. In real world, the star 
has to collapse, rearranging the angular momentum dis- 
tribution in its interior. What we can do, is determining 
the final stationary state (f), which by construction will 
be a stable rigidly rotating neutron star. We assume, that 
the transition conserves the baryon number of the star 
and is sufficiently rapid so that the angular momentum 
loss can be neglected. Therefore, the final stable configu- 
ration will have the same Mb and J as the unstable initial 
one, Mb,f = Mb,i , Jt = Ji. The difference between the 
parameters of these configurations (mass, equatorial ra- 
dius, moment of inertia) gives us the energy release, and 
changes in equatorial radius and rotation frequency, due 
to the collapse implied by the instability. 

Examples are shown in Figs llOli^ for our EOS model 
MUn. As the star becomes more compact, collapse is ac- 
companied by the decrease of the equatorial radius (by a 
few kilometers) and by a significant spin up. For a given 
EOS, the changes in radius, energy release, and spin-up 
are function of the angular momentum at the instability 
point: Ai?eq(Ji), A£'(Ji), A/(Ji). As wc see in Fig. El 
the energy release depends rather weakly on the rotation 
of the unstable configuration (i.e., on the value of Ji). It 
should be mentioned that our MUn model is only an ex- 
ample of the EOS resulting in the instability region within 
the hydrostatic equilibria. For EOSs with a weaker phase 
transition this instability region would be narrower and 
the changes of stellar parameters in the collapse would be 
smaller. However an approximate constancy of the energy 
release (i.e., its very weak dependence on J = Jj = Jf) 
seems to be a generic property of rotating neutron stars 
undergoing a collapse due to a first order phase transition. 



In order to discuss in more detail the energy release 
during collapses i — > f, we plotted in Fig. E| the gravita- 
tional mass of the star, M, as a function of angular mo- 
mentum, J, at fixed baryon mass: M = M{J)mj^ - Consider 
an initial configuration Ci. As the star looses angular mo- 
mentum, it moves down along line a, and reaches eventu- 
ally the cusp Ci (corresponding to the value of J = Ji). 
To continue moving on the dotted segment CiC2, the star 
would have to gain angular momentum and energy! 

As wc already mentioned, the evolution of the star 
beyond the instability point cannot be described by our 
hydrostationary model. The star can only collapse to the 
final configuration Cf , with the same values of Mb and J, 
i.e. along vertical arrow in Fig. 1121 Then, it evolves down 
the line b. 

We notice a very special role played by the point Cx 
at which line a and b cross. This is a degeneracy point, 
which corresponds to two very different configurations of 
the same Mb, M, and J. However, transitions between 
these two configuration are prevented by the huge energy 
barrier. 

The existence of the sharp cusps at Ci and C2 on 
the M{J)Mi, track is a very stringent test of the precision 
of the numerical code: it means that mass and angular 
momentum extrema (for fixed baryon mass) are reached 
exactly at the same point. This property follow s from the 
gener al relativistic relation dM = fid J + 7dMb (|Bardeen I 
Il972j) . Here, fl = 27r/, and 7 is the "injection energy per 
unit mass". This relation has to be strictly fulfilled by 
the stationary configurations. A graph, analogous to Fig. 
1121 can be plotted in the M — Mb plane for configura- 
tions with fixed angular momentum J. Also in this case, 
the existence of sharp cusps proves the correctness of the 
numerical code. 



7. Phase transition and pulsar timing and age 

It has been already pointed out bv ISpvrou fc StergioulasI 
(|2002ll , that the back bending phenomenon, resulting from 
the growth of a dense-phase core, can lead to significant 
difference between the actual pulsar age, ri, and that in- 
ferred from the measurements of the period P and pe- 
riod derivative, P, and denoted T2. The calculation of T2 
is based on quite strong assumptions. Firstly, pulsar ki- 
netic energy loss due to radiation is given by the magnetic 
dipole formula. Secondly, non-relativistic approximation 
is used, with pulsar kinetic energy given by ■^Ifi^, where 
il = 2'k/P and pulsar moment of inertia is constant, in- 
dependent of 

Following ISpvrou fc StergioulasI l)2002(l . we will use 
general relativistic notion of total pulsar energy, Mc^. 
Then the pulsar energy balance is 



dM 
"dT 



(2) 



where the right-hand-side is the pulsar energy loss rate via 
radiation of electromagnetic waves and particles. 
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Fig. 13. The relative mass-energy increase due to rota- 
tion of the star at fixed baryon mass, e(/) = — 
M(0))/A/(0)]A/b, for the EOS with a phase transition 
(MM model), is shown using solid line. Three solid lines 
are labeled by the gravitational mass of the non-rotating 
configuration (in solar masses). Dotted lines correspond 
to the e{f) curves calculated for the standard model, Eq. 
(|2Jl, with a = 4. 

Let us consider the increase of stellar energy due to a 
spin up to frequency / at constant Mb . In general relativ- 
ity, the increase is given by AM{f) = [M{f) - M(0)]A/b- 
In the standard model, we neglect the effect of / on stel- 
lar structure, so that AM = ij(27r/)^. This is a good 
approximation when rotation is slow and EOS is smooth 
(no phase transition) . Equation Q can be then rewritten 
as 

n = -^n'^-\ (3) 

where Iq = /(O) = const. 

In the case when angular momentum loss leads to a 
phase transition at the stellar center, the situation is much 
more complicated, because of the strong /-dependence of 
the pulsar structure in the vicinity of the phase transi- 
tion. This difference is illustrated in Fig. E| where we 
plotted the quantity e(/) = AM(/)/M(0) resulting from 
our calculations, and compared it with results given by 
standard non-relativistic model with constant I. As it can 
be seen in Fig. 1131 the approximation given by Eq. (|3J) 
significantly overestimates the change of frequency associ- 
ated with a given energy loss. This in turn can lead to an 
underestimation of the age of the pulsar - an example is 
presented in Fig. There, we plot the pulsar period, P, 
as a function of time, obtained by the integration of Eq. 
||5J) for magnetic dipole braking (a — 4). We also show 
results inferred from the observation of a 10 ms pulsar, 
assuming four selected values of the pulsar baryonic mass. 
Taking P = 10 ms and corresponding P the extrapolation 
backward in time using standard model, Eq. (|3Jl, diverges 
from exact results as soon as P < 5 ms. If the pulsar is 




-8 -6 -4 -2 



time [kyr/zCgg] 

Fig. 14. The evolution of the pulsar period when the en- 
ergy loss is described by the magnetic dipole braking with 
a = 4. Solid curves - results for our model MM with dif- 
ferent values of baryon mass. Dotted lines correspond to 
the standard model, Eq. The unit of time (horizontal 
axis) is 1000yrs/K28 where — k/10^* [cgs] is parame- 
ter entering Eq. ^ 

born with period of 2 ms, then its real age is significantly 
longer than P/2P. 

8. Discussion and conclusions 

The effect of phase transitions in dense neutron-star cores 
on the spin-down evolution of isolated stars was studied 
using analytical representations of EOSs. This guaranteed 
a very high precision of our 2-D calculations and allowed 
us to study a very large parameter space as far as the 
phase transitions themselves are concerned. In this way, 
we studied a very broad class of constant-pressure phase 
transitions associated with a density jump in the EOS. We 
also studied the case of transitions through a mixed-phase 
state. 

We limited ourselves to the case of normal rotating 
configurations, which are connected with the non-rotating 
ones by the angular momentum loss. We considered two 
types of instabilities which bound the sets of stationary 
configurations: mass shedding and instability with respect 
to the axisymmetric perturbations. The EOSs split into 
two sets: those producing a single family of stable station- 
ary configurations (static and rotating) of neutron stars 
and those producing two disjoint families of stable station- 
ary configurations. Conjectures concerning normal config- 
urations based on EOSs with a phase transition has been 
formulated. If an EOS yields two disjoint families of static 
configurations containing "twin neutron stars" of the same 
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baryon mass but different radius, then also stable nor- 
mal rotating configurations form two disjoint regions in 
the mass - equatorial radius plane, which contain "twin 
neutron stars" of the same baryon mass but of different 
compactness. 

Very often, neutron stars are called a second fam- 
ily of compact stars, the first family being composed of 
white dwarfs. Therefore, our conjecture means that an 
EOS which produces a third (disjoint) family of static 
compact stars, produces also a third (disjoint) family of 
rotating normal stars. 

We have also shown the existence of a very spe- 
cial class of "fine tuned EOSs" with phase transitions 
which produce marginally stable stationary configurations 
of normal neutron stars, which form a boundary separat- 
ing stable stationary configurations (a line in the mass- 
equatorial radius plane). 

Conditions on the density jump in constant-pressure 
phase transitions were derived, under which their presence 
in the EOS produces the back bending phenomenon in the 
spin down evolution. 

The case when a spinning down normal neutron star 
reaches an unstable configuration was studied in detail. 
The instability leads to neutron star collapse, associated 
with an energy release in a "corequake" , decrease of ra- 
dius, increase of central density, and spin up of the star. 
We have shown that the energy release associated with 
such a "corequake" depends rather weakly on the initial 
rotation frequency at the instability point. In our exam- 
ples, energy release was of the order of a few times 10^" erg. 

In the present paper we put accent on the numeri- 
cal precision and mathematical strictness. We hope that 
in this way we prepared ground for further studies of the 
impact of the phase transitions in dense matter on the 
structure, evolution, and dynamics of rotating neutron 
stars. These further studies will be performed using re- 
alistic EOSs available in the literature and taking into 
account important microscopic aspects of the phase tran- 
sitions. The kinetics of the phase transition coupled with 
stellar spin-down, and the ensuing neutron star corequake 
are now being studied. These topics will be the subject of 
our subsequent papers. 
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Appendix A: Analytical equations of state 

The relativistic polytrope, following iToooeil 1^^^, relates 
pressure p to baryon number density by 



(A.l) 



where F is the adiabatic index. The second coefficient, K, is 
often called the pressure coefficient ^. 



^ Unless otherwise mentioned, the coefficient K will be mea- 
sured in pc? /ff units, where p := 1.66 x 10^* g/cw?. and 
n ;= 0.1 fm-^ 



Dense matter is strongly degenerate, so that the T = 
approximation is valid. First Law of Thermodynamics implies 
then expression for energy per baryon 



7 - 1 



(A.2) 



where the constant C is set, for Wb = 0, to be the unit baryon 
rest energy^ C = niQC? . 

The energy density £ is thus given by 



K 



-nl + moc Tib 



F - 1 

The baryon chemical potential p is therefore 

, . p + £ TK r-i , 2 



Tib 



F - 1 



(A.3) 



(A.4) 



A.l. Constant pressure phase transition 

We assume that the phase transition stakes place in thermo- 
dynamic equilibrium. In a simplest case the transition from 
less dense pure phase A to high-density pure phase B occurs 
at constant pressure. Example of such transition is shown on 
Fie. IA.il Pure phases, A and B, will be approximated by poly- 
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Fig. A.l. Examples of phase transitions considered in the 
text; constant pressure phase transition (left), and the 
phase transition through the mixed-phase state (right). 



tropes, with adiabatic indices Fa and Fb, respectively. Let us 
also fix the mean baryon mass mj,^ and coefficient JsTa, and 
demand that the phase transition occurs between the baryon 
densities riA and tib, at constant pressure Pab and chemical 
potential p. From constancy of pressure and chemical potential 
we get 



^oA + Pab 



Fa 



nA(FA-l) nB(FB-l) 



(A.5) 
(A.6) 



A.2. Transition througii the mixed-phase state 

Lower-density phase A extends up to pressure Pi and density 
nh ~ ni, then follows mixed m-phase with volume fraction of 
dense phase xb increasing monotonically from at ni to 1 at 
nb = 712 and p = P2. Again, for simplicity, it will be assumed 



^ Following many authors, we put unit baryon mass equal to 
mass per nucleon in the ground state of atomic matter at zero 
pressure, which is ^^Fe crystal, mo = 1.66 x 10~^^ g 
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that the phase A and mixed-phase m can be approximated by 
polytropes (not a bad approximation, seeJBciEcr ct al. 2005). 
The parameters K^, Fm, and m^^^ (the mean particle mass in 
the mixed-phase) must be related to those of the A-phase in 
such a way that pressure and baryon chemical potential stay 
continuous across the A — > m transition point at n\, 



K 



A ■ ni 



Pi 



(Fa - l)(r,^ - 1) 



(A.7) 
(A.8) 



We will further assume that the high-density phase B, existing 
at nb > n2, p > P2, and p > P2, is pure quark matter, with 



at Ub > n2, p > t^2, and p > P2, is 
MIT bag model EOS dZdunikl I2OO0I1 . 

pi£) = ^{£-£o), n^ip) = mil +p/£f* 



(A.9) 



where £0 = poc^ and no are the mass-energy density and 
the baryon density of the quark matter at zero pressure. The 
baryon chemical potential of the B-phase is then equal to 



m(p) =Mo(l + 4p/£:) 



1/4 



(A.IO) 



where ^0 = £o/no. 

From the continuity of pressure, baryon density, and en- 
ergy density at the m — > B transition point, we get 

£0 = £2 - 3P2(n2) , 

no = n2/(l + 4P2/f2)'/'' . (A.ll) 
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